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The study of spin-related phenomena in materials requires knowledge of the precise form of effec- 
tive spin-orbit coupling for conducting carriers in solid- states systems. We demonstrate theoretically 
that curvature induced by corrugations or periodic ripples in single-layer graphenes generates two 
types of effective spin-orbit couplings. In addition to the spin-orbit coupling reported previously 
that couples with sublattice pseudospin and corresponds to the Rashba-type spin-orbit coupling in 
a corrugated single-layer graphene, there is an additional spin-orbit coupling that does not cou- 
ple with the pseudospin, which can not be obtained from the extension of the curvature-induced 
spin-orbit coupling of carbon nanotubes. Via numerical calculation we show that both types of 
the curvature- induced spin-orbit coupling make the same order of contribution to spin relaxation in 
chemically clean single-layer graphene with nanoscale corrugation. The spin relaxation dependence 
on the corrugation roughness is also studied. 



I. INTRODUCTION 

Graphene has attracted much interest due to its un- 
usual linear energy spectrum and electronic properties P 
It is also a promising material for spin-based applica- 
tions since spin relaxation (SR) in graphenes is expected 
to be weak due to the suppression of the hyperfine in- 
teraction of 12 C and the weak atomic spin-orbit coupling 
(SOC). Various spintronic devices^ and spin qubits^ on 
graphenes have been suggested. 

Recent experiments, however, measured SR times 
(SRTs) of the order of 0.1 — 0.5 ns, 4 which were much 
shorter than expected and comparable to other materi- 
als such as aluminuirP and copper. 6 Those experiments 
stirred up intensive theoretical investigation of the SR 
in single-layer graphene (SG) to clarify the origin of 
the short SRT. SG samples in RefW have low mobility 
~2000cm 2 /Vs and the SRT was found to increase with 
increasing carrier density, which is consistent with the 
Elliot- Yafet (EY) SR mechanism 7 arising from the inter- 
play of SOC and momentum scattering. However, it is 
not clear what kind of SOC is mainly responsible for the 
SR in SG, and various kinds of SOCs are examined. It 
was suggestecPthat ad atoms on SG can enhance the local 
SOC to the order of lOmeV, which results in short SRT 
comparable to the experimental value 4 for low mobility 
SG samples. Thus the adatom-induced SOC provides a 
possible origirPE of the short SRT in Refffl If this SOC 
is indeed the main origin, which requires further confir- 
mation, it implies that the SRT may be enhanced in SGs 
without adatoms. To find out how greatly the SRT can 
be enhanced, various factors limiting the SRT should be 
examined carefully. In exfoliated SGs, charged impuri- 
ties or surface phonons in the substrate^ can induce the 
effective SOC of the order of 0.01 meV. Also, effective 
SOC of the order of 0.01 - 0.1 meV^ can be induced 
by local curvature effects arising from corrugations^^ 



and periodic ri pples^ ^ observed in suspende d 12 * 17 * 18 * 
and exfoliated ^ 16 * 19 * SGs. Each substrate-induced and 
curvature- induced S OC co mbined with momentum scat- 
tering off impuritie d 11 * 20 * or electron-electron Coulomb 
scatteringpl and substrate-induced or curvature-induced 
SOC itself 22 can cause SR. The estimated SRTs are at 
least of the order of 10 ns, which is about two order of 
magnitude larger than measured SRTs. 4 

For a more reliable investigation of spin-related phe- 
nomena such as SR and spin coherence, it is important 
to know the precise form of SOC. In this respect, existing 
studies of the curvature-induced SOC in SGs are not sat- 
isfactory, since the precise f orm o f this SOC is not known. 
Previous theoretical studies^^IMI] assumed a particular 
form of the SOC inferred from that of a carbon nan- 
otube (CNT)P However the SOC-induced energy split- 
ting measurement^ in ultr a- clea n CNTs revealed that 
the form of the SOC in Ref ! 1Q | 241 does not provide a sat- 
isfactory description of the SOC effect in a CNT. Later 
theoretical works^L^ reported that the correct form of 
the SOC in a CNT has an additional term in addition to 
the SOC form in Ref ! 1Q * 24 ( The additional term does not 
couple to the pseudospin degrees of freedom in a CNT 
and thus differ s qualitatively from the previously known 
SOC termP^H w j 1 j c ] 1 CO uples to the pseudospin. It was 
demonstrated in Ref.^that the interplay of the two SOC 
terms can explain the SOC-induced energy splitting mea- 
surement pM3 it i s then reasonable to expect that a sim- 
ilar additional SOC term may exist for SGs as well. 

Finding the correct curvature-induced SOC can have 
an implication not only on SR in the current experi- 
ments, 4 but also on SR in chemically clean SGs. As 
technology for the preparation of SGs progresses so as 
to diminish impurity effects, the fundamental limit of SR 
would be governed by intrinsic sources. Recalling that 
atomically flat SGs are thermodynamically unstable,^ 
the curvature-induced SOC is expected to be one of the 
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intrinsic origins determining the upperbound of the SRT 
in chemically clean SGsP^ We believe that this expec- 
tation is very plausible since SGs and CNTs are very 
similar chemically and th e dominant source of the SOC 
in ultra-clean CNTs^^Zlis the curvature-induced SOC. 

In this paper, we show that in addition to the 
curvature-induced Rashba-type SOC reported in an ex- 
isting theory,^ which couples to the pseudospin, there 
exists an additional type of the curvature-induced SOC, 
which does not couple to the pseudospin and thus ap- 
pears as the diagonal term in the pseudospin representa- 
tion. In this respect, this additional S OC term is similar 
to the additional SOC term in C NTs PEH However these 
two additional terms have topologically different origins 
since the periodicity in the circumferential direction aris- 
ing from the tube topology of CNTs is absent in SGs. 
Thus the additional SOC term in SGs cannot be obtained 
from the additional SOC term in CNTs through a trivial 
extension. This point was not addressed in an earlier at- 
tempt 26 to derive the curvature-induced SOC in SGs. In 
this paper, we also examine the SR due to the curvature 
in chemically clean SGs and show that the effect of the 
additional SOC on the SR is comparable in magnitude 
to the effect of the previously known curvature-induced 
SOC. Thus the additional curvature-induced SOC term 
should be included in studies of spin-related phenomena 
in SGs. We also investigate SR dependence on the frac- 
tal dimension of the corrugation roughness j 14 * 30 * and show 
that corrugation roughness affects the SR. 

This paper is organized as follows. In Sec. [TTJ we 
show analytical expressions of two types of the curvature- 
induced SOCs in SGs, and then demonstrate how to 
obtain them microscopically using tight-binding (TB) 
Hamiltonian of the local curvature effect and the atomic 



SOC of carbon. In Sec. Ill we investigate effects of the 
effective SOC on SR based on the kinetic equation of the 
spin density operator in corrugated SGs, and calculate 



numerically SRT. We give a brief summary in Sec. IV 



K 2 = — Ki, 7r-band states with wave vector k = (k x ,k y ) 
[|k| <C |K!( 2 )|] relative to the Ki( 2 ), can be described by 
the effective Dirac Hamiltonian Hoirac- When the two 
sublattices A and B of the honeycomb lattices are used 
as bases, the Hoirac is written as^H 



e iT * a (r z k x -ik y ) 
Y (T z k x +ik y ) 



(1) 

where vp is the Fermi velocity, and r z = ±1 denotes the 
K!( 2 ). Here, the angle a is defined counterclockwise from 
the y direction and the C-C bond vector b cc -direction 
(Fig.0). 

Local curvature effects arising from smooth corru- 
gations or ripples in SGs combined with atomic SOC 
of carborP^ generate the effective SOC for the 7r-band 
states. 10 The corrugated SG can be described as a sur- 
face z — h(r) with height h(r) as a function of spatial 



position r 



in the two-dime nsiona l (2D) refer- 



ence plane (Fig. 0). Experiments^"^^ suggest that 
h(r) shows fluctuations of the order of 1 nm over scales 
~10nm, which can be specified by hihj <C 1 (i, j = x,y) 
with hi being the partial derivative along the i direction. 
When the local structure of the SG surface has two prin- 
cipal curvatures fti(r) and ^(r), we find that the ^ soc ( r ) 
is written as 



^soc( r ) 



Cr z /jb(r) • s Ce lTza i/(r)(r z s y + is x ) 

(e- iTza is(r)(T z s y - is x ) C'TsMM ■ s 

(2) 



II. EFFECTIVE SPIN-ORBIT COUPLING 



where s = (s x ,s y ,s z ) denotes the real spin Pauli ma- 
trix. Here, £ and are material parameters, and /x(r) 
and v(y) are geometrical parameters containing the lo- 
cal curvature information. For a = 0, the off-diagonal 
part of Eq. (pi) reduces to the Rashba-type SOC oc 
(& y s x — T z a x Sy) reported in existing theorief^^, where 
a = (cr x ,cr y ,cr z ) denotes the sublattice-pseudospin Pauli 
matrix. On the other hand, the diagonal part of Eq. ([2| 
is one of main results of this paper, which was not ignored 



In this section, we present the effective Dirac Hamilto- 
nian in flat SGs, and then, the effective SOC Hamiltonian 
arising from the local curvature effect combined with the 
atomic SOC in corrugated SGs. After that, we demon- 
strate how to evaluate the effective SOC using the TB 
Hamiltonian via second-order perturbation theory, and 
why the SOC can not be inferred from the extension of 
that of CNTs. 



A. Analytic expression 




Graphene honeycomb lattice with two sublattices A 
and B (Fig. [I]), has 7r-band consisting of p z orbitals 
near the Fermi level situated at hexagonal corners in 
the Brillouin zone. Close to the Dirac point Ki and 



FIG. 1. (Color online) Single- layer flat graphene honeycomb 
lattice with two sublattices A and B. The angle a is defined 
counterclockwise from the y to the C-C bond vector b cc . a = 
I b cc | ~ 1.42 A is the carbon-carbon distance. 
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in previous studies. Note that this additional SOC term 
in the diagonal part does not couple to the pseudospin. 

The material parameters ( and ( f are, respectively, 
given by ( = a\ so (e p - e 8 )(V£ + V^)/{12V S % 2 ) and (> = 
aX soV^ p /(2(V^ p - V p a p )), where A so is the atomic SOC 
strength of the p orbitals, e s ^ is the atomic energy for 
s(p) orbitals and a is the carbon-carbon distance. Here, 
and Vpp^ represent the coupling strength in the ab- 
sence of the curvature for the a coupling between nearest- 
neighbor s and p orbitals and the ir(a) coupling between 
nearest-neighbor p orbitals, respectively. 

The geometrical parameters ^(r) and ii(r) are given 
as follows, ^(r) is a sum of the two local curvatures, 
^(r) = fti(r) + k,2(t) and /x(r) is a weighted sum of the 
two local curvatures, 



/x(r) = (fj, x (r),fj, y (r)) 

2 

= Kj(r) (sin [20-,- (r) - 3a] , cos [20,(1-) - 3a]) , 



(3) 



where 0i(r) and 02 (r) are the angles between a local 
x / (r)-direction and the two principal curvature directions 
(see Fig. [5] for their precise definitions). Since the two 
principal curvature directions are mutually orthogonal, 
2 (r) =0i(r)±7r/2. 




FIG. 2. (Color online) Schematic of a partial convex struc- 
ture in a corrugated SG surface described by z = h(r) where 
r= (x,y) within the 2D reference x?/-plane. The z'(r) is the 
unit vector normal to the local tangential plane and the h' cc (r) 
(green arrow) is the projected vector from the local C-C bond 
onto the tangent plane specified by the (blue) dashed quad- 
rangle normal to the z'(r). The y x (r)-direction in the tangent 
plane is defined to make the angle a [Fig. [I] clockwise from 
the bc C (r) in the tangent plane. Then the x / (r)-direction is 
determined automatically to be mutually orthogonal to z x (in- 
direction and y / (r)-direction. Dotted (blue) lines are cur- 
vature lines tangent to principal directions along orthogonal 
(blue) arrows where the direction of arrowheads is arbitrary. 
The angle 0i(2)(r) is defined by x / (r)-direction and the princi- 
pal directions, and 02 (r) = 0i(r) ±7r/2. The sublattice A and 
its nearest-neighbor sites Bi (I — 1,2,3) with two curvature 
lines are indicated. 



Both diagonal and off-diagonal terms in Eq. (2j) are 
invariant under the time reversal symmetry. The off- 
diagonal term is allowed when the mirror symmetry with 
respect to the x?/-plane is broken.^ On top of that, owing 
to the broken C3 symmetry in corrugated SGs, the diago- 
nal term of the form r z s x ^ is allowed as well. Note that 
the diagonal term has the 0i( 2 )(r)-dependence. However, 
it can not be obtained from the naive extension of the 
diagonal SOC of CNTs that depends on the chiral angle 
through cos(30)P^Sl Quantitatively, both terms have 
the same order of magnitude. Using A so « lOmeVp^ 

4.20 eV, V p a p = 5.38 eV, 



-7.3 eV, e p = OeV, V & 



y-Tr 
pp 



sp 



2.24 eV in RefP we obtain (' « 0.21meVnm 
and £ ~ 0.15meVnm. 

The 0j(r) and the a-dependence of H soc (r) combined 
with the SG morphology have interesting implications. 
First of all, since 62 (r) = 0i(r) ± 7r/2, the diagonal 
term depends on the difference between two principal 
curvatures in contrast to the off-diagonal term depend- 
ing on their sum. As a result, at umbilical points where 
fti(r) = /^(r), the diagonal term disappears while the 
off-diagonal term remains, which confirms that every tan- 
gent vector can be a principal direction at the umbilical 
points. 35 In contrast with the umbilical points, at sad- 
dle points where fti(r) = — ^(r), the diagonal term only 
appears while the off-diagonal term disappears. 

In addition, for periodic ripples having unidirectional 
curvature direction (fti(r) 7^ and ^(r) = 0), the diago- 
nal SOC field direction varies depending on the principal 
curvature direction and the honeycomb lattice orienta- 
tion as well. When a = q/3 (q G Z) so that the x is 
along the zigzag direction [Fig. fl], for example, the di- 
agonal term is ex Ki(r)s y for 0i(r) = and 7r/2, while 
it is oc Ki(r)s x for 0i(r) = 7r/4. On the other hand, 
when a = q/3 + 7r/6 so that the x is along the arm- 
chair direction [Fig. HI, the diagonal term is oc Ki(r)s x for 
0i (r) = and 7r/2, while it is ex Ki(r)s y for 0i(r) = 7r/4. 
Recent experiments shows that the periodic ripple line 
direction can be controlled^ and the multiple periodic 
ripple domains with different ripple lines can occur 
which implies that the variation of the diagonal term de- 
pending on the curvature direction and the honeycomb 
lattice orientation could be significant for the analysis of 
spin transport in SGs with those structural deformation. 
After all, in order to describe precisely the curvature- 
induced SOC in SGs both terms should be considered on 
an equal footing. 

The / H S oc( r ) can be written as an explicit function 
of h(r) in the SG surface, z = h(r) where hihj <C 1 
(i, j = x,y). The principal curvature is the eigenvalue of 
the shape operator S = T^ X T2 with ^1(2) being the first 
(second) fundamental form.^ Moreover, since its eigen- 
vector #1(2) (r) is along the principal direction, 0i(2)(r) 
can be determined by the V\{2) ( r ) an d x'(r), which is as- 
sumed as x'(r) w (1,0, h x ) in the zeroth order of hihj. 
T\ and T2 are, respectively, given b}E^ 
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J 7 !- 



1 + hl 
h x hy 1 



h x hy 



, ^2 



(4) 

and then, the fti( 2 )(r) and #i( 2 )(r) are written as 

*i(2) (r) « (/^ + ^ - (-1) 1(2) - ^) 2 + 4hly] * ) /2 
andvi( 2 )(r) w [ai( 2 )X+^y+(ai( 2 )ft a; +/8ft 3/ )z]/jV, respec- 
tively, where a 1(2) w ft^ - h yy - (-1) 1(2) [(^ - ft^) 2 + 
4ft 2 ] 2 and /3 « 2ft X2y with normalization constant A/". 



B. Tight-binding Hamiltonian 

Now we demonstrate how to derive Eq. ([2| microscopi- 
cally. Using second-order perturbation theorjJ*^ with the 
local curvature and the atomic SOC as weak perturba- 
tion, we obtain Eq. (J2|. Thus we need to know the TB 
Hamiltonian of the local curvature H c and the atomic 
SOC H so taking the honeycomb lattice structure into ac- 
count exactly. Most parts of this subsection are devoted 
to the TB theory of the ir-a hybridization among s and 
p orbitals. The unoccupied d orbital effects^^ are dis- 
cussed at the end of this subsection. 

Firstly, H c is determined purely by local curvature ef- 
fect dependent on fti( 2 )(r) and #i( 2 )(r) at three bonds be- 
tween an atomic site r A in the sublattice A and nearest- 
neighbor sites B\ (I = 1, 2,3) in the sublattice B [Fig. [2]. 
For the smoothly corrugated SGs, we consider solely one 
of two principal curvature effects assuming that the other 
principal curvature effect is absent. There are also mu- 
tual curvature effects between them, but those effects 
can be ignored for the smooth corrugated surface with 




"«i( r >0 




(b) 



FIG. 3. (Color online) Schematic of a partial convex SG sur- 
face with sublattices A and B\ (I = 1, 2, 3) indicated in Fig. [2] 
which is assumed to have unidirectional principal curvature. 
Although the original surface [Fig. [2] has finite principal cur- 
vatures ki(ta) and ^(^a), each curvature effect is evaluated 
under the assumption that the other principal curvature is 
zero. ki(ta) / and ^(fa) = in (a) and ki(t a ) — and 
^2(1* a) 7^ in (b). Dotted (blue) lines represent curvature 
lines, and principal direction at is tangent to the curva- 
ture lines. The arc length (j — 1,2) between a solid (green) 
line passing A atom normal to the principal direction, and its 
parallel solid (green) line passing B\ atom, is related with uj\ 



by ^«^(r)/2 [Eq. §], 



hihj <C 1 = x,y) since they are higher order in 

a/lZ 1(< 2)(r A ) where Ui(2)(r A ) = k~^(t a ). Then, it is 
enough to deal with two unidirectional ripple structures 
whose principal directions are orthogonal as shown in 
Fig. [3j For the unidirectional ripple structure, each prin- 
cipal curvature and the angle between the principal direc- 
tion and the x'(r) at three neighbor C-C bonds around 
the i\4 are same as ^i( 2 )(r^) and 0i(2)( t a), respectively. 
Then, up to first order in a/7^i( 2 )(r^), the H c can be 
approximated as H c « Hi + H2 , each of which describes 
one principal curvature effect in the absence of the other 
one. Two principal curvatures yield correlation effects 
between term in H C1 but they are higher order in term of 
a/7^ 1 ( 2 )(r) and thus can be ignored. The Hj (j = 1, 2) is 
written a^P^H 



= [- 



YA 1 = 1 77=t,4- 



a f x' 

1 \^V A 



t r z' 
iV Birj 



1 j y-r A V Bv 



TAV BlT] J 



V Birj 



^VAV BlT] , 



U YAT)^B { 



H.c, 



(5) 



where cZ 



and 



represent the annihila- 



tion operators for p x > , p y > , p z > orbital states with eigen- 
spinor Xt(l) a l° n g z' a ^ a carbon atom r. Here, Sp 
Xj l , and Yj l represent the curvature-induced coupling 
strength of p x r, p y > , and s orbitals with a nearest- 
neighbor p z t orbitals. Their expressions are written as 



■AxlcaaiOjiTA)] - JY\ sin[^(r A ) 



and y/ f =o;.Xjsin[^(rA)]+w -y/cos[^(rA)], with 



a/(2n j (r A ))sm[0 J (r A )-a] 
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a/(2^(r A ))sin[7r/3- 



Oj (r A ) + a] , and « a/(2Kj (r A )) sin[7r/3 + Oj (r A ) - a] . 
Here, S\, X\, and Yj are, respectively, written as 



Sj = VZ,sm\fi j (r A )], 



3 

S? = V£cos 



j - ~V; P - V P 1 cos' \Pj(t a )] - Vp p sm 2 \fij(r A )] , 



x] = -v; p -v; pCOS 2 



X s = V n + V v cos 2 
3 pp K pp w ° 



7T 

L6 



-hi* a) 

■Pi<? a) 



7T 

3 

7T 

L6 " 



hi* a) 



Yf = j (V p l-V; p ) sin [2ft- (r A )], 



^ 2 = .(^-^,)sin 



y 3 = - (v n - V° 

3 O VPP PI 



2Pj ^a) 



2tt 



(6) 



where (3j(r A ) = 6j(r A ) — a. Note that the local curvature 
generates interband transition between 7r-band and a- 
band consisting of s, p f x , and p' y orbitalsF 10 * 26 * 39 * 
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Secondly, H so is given by H so = A so ^ r L r • S r pSESl 
where L r (S r ) is atomic-orbital (spin) angular momentum 
of an electron at r, and can be expressed b} * 10 * 26 * 

11 so — ^ \ ri C rt C rt C r± T" ~h fcC^ 

+<; t 4-<; t <)+H.c. (?) 

Note that on-site coupling between p z > and jv^/) orbitals 
generates ir-cr interband transition that, combined with 
the curvature- induced interband transition, contributes 
to the H S oc(r). 

In CNTs, however, the coupling between p z > orbital 
and p orbital along the CNT axis has no contribution 
to the effective SOC. For instance, if the Eq. ^ is ex- 
pressed with eigenspinor along the CNT axis parallel to 
y^x^-direction instead of XtU)> the coupling between 
p z > and p y '(x') orbitals acquires a phase term e ±% ^ with 
the azimuthal angle (j) along the circumference^^, and, 
as a result, its contribution averages out to zero after 
integration over the circumference Hence, the trivial 
extension from the curvature-induced SOC of the CNT 
to that of the SG is inevitably incomplete because of their 
topological difference in the geometrical structure. 

Second-order process resulting from two consecutive 
7T-(T interband transitions gives rise to the effective 
SOC Hamiltonian for the unperturbed 7r-band states of 
Eq. QP^For the cr-band calculation, we use the Slater- 
Koster parametrization 42 of s, p xi and p y orbitals. Dur- 
ing the second-order process, pseudospin-conserving and 
pseudospin-flipping processes occur; the former and the 
latter generate, respectively, the diagonal and the off- 
diagonal SOC in H S oc(r) [Eq. QjPThe resulting effec- 
tive SOC, then, is expressed in local coordinate axes of 
x'(r), y'(r), z'(r). However, for the smooth corrugation 
where hihj <C 1 (i, j = x : y), x'(r), y'(r), z'(r) in the re- 
sulting effective SOC may be replaced by x, y, z to obtain 

Hsoc( r ) hi Eq. pf since the difference between the two 

— _ i 

sets of axes [for instance, z/(r)-z = (l + |V/i(r)| 2 ) 2 w 1] 

generates higher order terms in a/7^ 1 ( 2 )(r). 

Lastly, we remark d orbitals effects on the curvature- 
induced SOC, which were not considered in our calcula- 
tion. For the spin-orbit induced gap, the d orbital effects 
are recently reported to be important since the hybridiza- 
tion between d orbitals and p orbitals generate the gap, 
which is linear in the atomic SOC strength A^ of d or- 
bitals.^^This linear contribution can be dominant over 
the p orbital contribution, which is quadratic in A so . In 
contrast, the contribution of the d orbitals to the exter- 
nal field-induced Rashba SOC is smaller than that of the 
p orbitals. 38 The coupling between p and d orbitals can 
also lead to the curvature-induced SOC in the first order 
of a/7^i( 2 )(r). Due to the symmetrical reason, however, 
its form is expected to be same as the form of H soc (r) 
[Eq. ([2])] except parameters in ( and (' . The curvature- 
induced SOC owing to the d orbitals can be generated as 



follows. The 7r-band consisting of the p z orbital and the 
nearest-neighbor d xz and d yz orbitals can couple with the 
<j bands through on-site coupling between d orbitals in 
7r-band and cr-band by the atomic SOC. Then, electrons 
in cr-band can delocalize to the nearest-neighbor p and d 
orbitals within a bands. Finally, the p and d orbitals in 
cr-bands can couple with the nearest-neighbor p z orbitals 
due to the curvature, yielding the curvature-induced SOC 
that is proportional to a/7^i(2)(r), A^, and (e^ — e p )~ 2 . 
Based on the estimation of the Rashba SOC strength 
and comparison with the p orbital contribution given in 
Ref! 3 -^, we can compare roughly the d orbital contribution 
of the curvature-induced SOC with the H soc (r) [Eq. [2]. 
The Xd is smaller than A so , and further, (e^ — e p )~ 2 mul- 
tiplied by hopping parameters associated with d and p 
orbitals is a small quantity compared with parameters in 
same dimension in ( and (' . Hence, we believe that the 
contribution of d orbitals to the curvature-induced SOC 
is smaller than that of p orbitals. In order to get more 
reliable quantitative estimates, systematic study consid- 
ering coupling between p, s orbitals and d orbitals is re- 
quired. 38 



III. SPIN RELAXATION 

In this section, we investigate effects of the geometric 
curvature on SR in chemically clean SGs with nanoscale 
corrugations. Motivated by recent experiment^ about 
the corrugation roughness of SGs, we study SR depen- 
dence on the fractal dimension of the corrugation rough- 
ness within the assumption that both SOC and mo- 
mentum scattering 30 arise mainly from the corrugation. 
Since the substrate-induced SOC effect on SR was ad- 
dressed in a previous theory, 9 we focus on the curvature- 
induced SOC here. 



A. Kinetic equation of the real spin density 
opeator 

In order to calculate SRT, we use the kinet ic equa- 
tion of the real spin density operator p ! 22 * 43 1 Besides 
the % SO c( r ) 7 we include the strain-induced vector po- 
tential H v (r) that occurs in corrugated SGsJ^HH Al- 
though H v (r) is independent of the spin, it is still rel- 
evant for the SRT since it affects the momentum scat- 
tering rate. Then the local potential Hamiltonian H p (r) 
is written as H p (r) = H soc ( r ) + %v( r ) which is consid- 
ered as weak perturbation in the kinetic equation. We 
assume (H p (r)) w where (• • • ) denotes ensemble aver- 
age. Since H p (r) is expected to change slowly over scales 
larger than the lattice spacing so that the hybridization 
between r z — 1 and r z = —l states can be ignored, we cal- 
culate SRT for the electron state of the Hoirac [Eq. ([!])] 
when r z = l and a = 0, where eigenenergy and eigen- 
state I -0k) are > respectively, given by = hvpk with 
k = |k|, and |^ k ) = l/V^e - ^ l) T |k) with polar angle 
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<£>k of the wavevector kP^ Then a 2 x 2 real spin density 
operator p^ diagonal in momentum has a representation 
in eigenspinors chosen to be along z direction. The ki- 
netic equation of pk is written a d * 43 * 

^ + J[Wk]=C k , (8) 

where Vkk' = (^kl^pO*)!^'), and Ck is the collision inte- 
gral describing momentum scattering, which is given by 

Ck = 7r/^Ek' [2Vkk'Pk'Vk'k - {pk, Vkk'Vk'k}] S(e^ - £k). 
The commutator in Eq. Q can be ignored because H p (r) 
has no regular contribution. Since the energy scale of 
T-L v (r) is typically larger than that of H soc (r)^ the SRT 
would be longer than a time scale in which the momen- 
tum distribution during momentum scattering becomes 
isotropic.^ Thus, the pk can be represented as Pk = P+/4 
with the anisotropic part due to the ?{ soc (r), where the 
bar means averaging over cp^. Inserting the Pk = P + Pk 



and l/r| = 1/rJ + 1/r*, where J and J' are, respec- 
tively, defined as J = J( 2 (kR) + Jf(kR), and J' = J$(kR) 
with Jq(i)(x) being the Bessel function of zeroth (first) 
order and R = |r — r'|. If we assume (p x (r)p x (r f )) « 
(p y (r)p y (r f )) and (p x (r)p y (r f )) « for the random cor- 
rugation, we obtain « and « 2r|, which 
means that spins out of plane relax twice as fast as spins 
in plane. 20 Note that since the carrier density n is given 
by n = k 2 /V, the SR vanishes at the Dirac points where 
n = OP 



Figure |] shows SRTs, (t§) that are calculated nu- 
merically in corrugated SG surfaces for a = and are 
ensemble- averaged for ten surfaces. Those surfaces are 
approximated as square lattices, which can be justi- 
fied since the H soc ( r ) contains the curvature effects tak- 
ing into account the honeycomb lattice. The surface 
with fluctuating h(r) is constructed based on the height- 
height correlation function g(R) = ((h(r) — h(r')) 2 } — 
2j 2 (l — e~( R /^ ) , where H characterizes the fractal di- 
mensi onj 14 * 30 * 49 ! £ represents the correlation length, and 
7= y/ (h 2 (r)). As the fractal dimension H increases, the 
corrugated surface becomes smooth. 

A recent experiment!^ in suspended SGs shows that 
the height of corrugations grows approximately linearly 
with increasing length, which corresponds to a corru- 
gated surface for 2H = 2. For exfoliated SGs on the 
substrate, on the other hand, the fractal dimension with 



into Eq. ([8]), we obtain the equation of p as, 

^^jd^jd^lv^w^ - \ [p, K k ,\ 2 f) , (9) 

where d = k/^hvp) is the density of states. Note that 
collision term consisting of the product of H soc ( r ) and 
Hv(r) disappears since there is no correlation between 
them. More specifically, since the vector potential is writ- 
ten as quadratic in hij j 44 l 45 l the collision terms propor- 
tional to (hjhj/j,(r)) or (hihji/(r)) vanish by the Wick 
theoremPSl Thus, the random gauge field does not con- 
tribute to the SR in chemically clean SGsP^ 

Since the real spin density P = (P x , P yi P z ) can be eval- 
uated by P = Tr (ps), we can obtain the kinetic equation 
of P by tracing Eq. ([9| multiplied by s over the real spin 
space, which is written as dPi/dt ~ —Pi/rf (i = x,y,z) 
with rf being the SRT of spin along the z-direction. Here, 
^~x(y) an d r z are evaluated as 



(10) 

I 

2H « 1 was observed.^ There is a controversy over the 
origin of observed corrugation roughness!^ Here, we in- 
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FIG. 4. (Color online) Spin relaxation time (r|) as a function 
of the carrier density n — k 2 /it and the fractal dimension H. 
The (t^ ) are ensemble-averaged quantities over ten surfaces 
for a = 0. The corrugated square lattice with the height h(r), 
the lattice constant of 0.2 nm, and the side length of 200 nm, 
are constructed numerically based on the random midpoint 
displacement methocP^ with height-height correlation func- 
tion g(R) = 27 2 (l - e~ w ° 2H ). H ere w e set 7 = 0.3nm and 
£ = 20 nm from recent experiments^^ and interpolate dis- 
crete h(r) in order to mimic the continuous SG surface, (a) 
(Tz) as a function of n for 2H — 1 (red solid line) and 2H — 2 
(blue dashed line), (b) (t|) as a function of n for 2H — 1 
considering ?4oc(r) (red solid line) and only the off-diagonal 
SOC in 'Hsoc(r) (black dashed line), respectively. 
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vestigate SR in corrugated SGs for 2H = 1 and 2H = 2 
representatively. As shown in Fig. [4] (a), the (r|) becomes 
shorter with decreasing 2H, which can be understood in 
terms of a growing irregularity of H soc ( r ) [Eq.|2] induced 
by rougher surface corrugation. Also, we check that the 
( T x{y)) * s approximately twice as large as the (r|) (not 
shown here). 

Next, in order to investigate the effect of the diagonal 
SOC in % SO c( r ) 7 we compare two SRTs, one calculated 
considering the entire H soc (r) and the other considering 
only the off-diagonal SOC term in H soc ( r )? respectively, 
in SGs with 2H = 1. As shown in Fig. |] (b), the lat- 
ter SRT is noticeably larger than the former true SRT. 
Hence, it is necessary to consider both SOC terms for 
the precise analysis of SR depending on the corrugation 
roughness. 

The calculated SRT as shown in Fig. [4] is at least of the 
order ~ 10 ns, and becomes shorter monotonically as the 
carrier density n increases, which is primarily because the 
density of states d increases linearly with n 1 / 2 [Eq. (10)]. 
We check that the integrand of Eq. ( 10 ) has little depen- 



dence relatively upon k. This implies that the effective 
spatial range of the random SOC affecting the SR sub- 
stantially does not depend significantly on the variation 
of k over the scale ~ lOnrn -1 . Our quantitative and 
qualitative results are in agreement with a recent the- 
ory.^ As mentioned in the Ref! 2 ^, however, those results 
cannot explain recent spin transport experiments.^ 

One possibility for that discrepancy is that the ef- 
fect of charged impurity in the substrate causing mo- 
mentum scattering 51 is not considered in our calculation. 
The momentum scattering combined with the curvature- 
induced Rashba-type SOC was already investigated the- 
oretically,^ but it cannot explain the measured SRT, ei- 
ther. However, the inclusion of the additional diagonal 
SOC in H S oc( r ) hi addition to the off-diagonal SOC into 
the EY mechanism could give a possibility to calculate 
the SRT more precisely. In the EY mechanism^ spin flip 
scattering occurs at the very moment when the momen- 
tum scattering takes place since the electron wave func- 
tions normally have an admixture of the opposite-spin 
states due to the SOC. 



IV. CONCLUSION 

We have demonstrated that the combined action of the 
curvature and the atomic SOC of carbon gives rise to two 
types of the effective SOC in SGs with corrugation or pe- 
riodic ripples. One of them couples with the sublattice 
pseudospin, which corresponds to the Rashba-t ype SO C 
in a corrugated SG reported in previous theories } 10 I 2 °P^1. 
The other SOC, on the other hand, does not couple with 
the pseudospin, and was not recognized in previous litera- 
ture. The additional curvature-induced SOC depends on 
the principal curvature direction, which is similar to the 
curvature-induced SOC of CNTs whose diagonal term 
in pseudospin representation depends on the chiral an- 
gle of CNTs. However, the curvature-induced SOC in 
SGs can not be obtained from the trivial extension of 
the curvature-induced SOC in CNTs due to their dis- 
tinct topological structure between the SG surface and 
the CNT cylinder. 

We have also investigated SR in chemically clean SGs 
with nanoscale corrugation, and found that the diagonal 
SOC makes the same order of contribution to SRT as 
the off-diagonal SOC. The SRT becomes longer as the 
fractal dimension of the corrugation roughness increases 
since the irregularity of the SOC decreases in smoother 
SGs. The calculated SRT in chemically clean SGs, how- 
ever, can not explain recent experimental results of SR 
in current exfoliated SG samples. 

A natural direction for future research would be to 
calculate SRT in the presence of charged impurities that 
cause momentum scattering} 5 -^ considering both the diag- 
onal and the off-diagonal SOC. In addition to the SR in 
SGs, we expect that the curvature-induced SOC [H soc ( r )] 
may be relevant for the analysis of other spin-related phe- 
nomena in SGs. 
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